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Abstract. The Gerasimov-Drell-Hearn (GDH) sum rule connects the anomalous contribution to the mag- 
netic moment of the target nucleus with an energy-weighted integral of the difference of the helicity- 
dependent photoabsorption cross sections. Originally conceived for real photons, the GDH integral can be 
generalised to the case of photons with virtuality Q 2 . For spin-1/2 targets such as the nucleon, it then 
represents the non-perturbative limit of the first moment A of the spin structure function gi(x,Q 2 ) in 
deep inelastic scattering (DIS). The data collected by HERMES with a deuterium target are presented 
together with a re-analysis of previous measurements on the proton. This provides an unprecedented and 
complete measurement of the generalised GDH integral for photon-virtuality ranging over 1.2 < Q 2 < 12.0 
GeV 2 and for photon-nucleon invariant mass squared W 2 ranging over 1 < W 2 < 45 GeV 2 , thus covering 
simultaneously the nucleon-resonance and the deep inelastic scattering regions. These data allow the study 
of the Q 2 -dependence of the full GDH integral, which is sensitive to both the Q 2 -evolution of the reso- 
nance form factors and contributions of higher twist. The contribution of the nucleon-resonance region is 
seen to decrease rapidly with increasing Q 2 . The DIS contribution is sizeable over the full measured range, 
even down to the lowest measured Q 2 . As expected, at higher Q 2 the data are found to be in agreement 
with previous measurements of the first moment of gi . From data on the deuteron and proton, the GDH 
integral for the neutron has been derived and the proton-neutron difference evaluated. This difference is 
found to satisfy the fundamental Bjorken sum rule at Q 2 = 5 GeV 2 . 
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1 Introduction 

The Gerasimov-Drell-Hearn (GDH) sum rule connects an 
energy-weighted integral of the difference of the helicity- 
dependent real-photon absorption cross sections with the 
anomalous contribution k = ^^r 1 — Z to the magnetic 

• • [TTi 

moment /i of the target nucleus with atomic number Z M 
(or nucleon Jl];||]) : 

I [^M-^M]— = ~j^-« a . (i) 

Here cr^ and are the photoabsorption cross sections 
for relative orientation of the photon spin anti-parallel and 
parallel to the nucleus spin /, v is the photon energy in the 
target rest frame, v is the photoabsorption threshold, M t 
is the nucleus mass, a the electromagnetic fine-structure 
constant and e the elementary charge. This sum rule pro- 
vides an interesting link between the helicity-dependent 
dynamics and a static ground state property of the target 
nucleus. 

The GDH sum rule holds for any type of target, i.e. it 
is valid for protons, neutrons or nuclei. It is also consid- 
ered to be important in electroweak physics |l| . The GDH 
sum rule is derived starting from the Compton forward- 
scattering amplitude following the general physics princi- 
ples of Lorentz and gauge invariance and is non-pertur- 
bative in nature. The only questionable assumption in its 
derivation is the use of an unsubtracted dispersion rela- 
tion. For the proton (k p = +1.79) the GDH sum rule 



prediction is — 204/ib, for the neutron (k„ = —1.91) it is 
—233 /ib. The prediction for the deuteron (k^ = —0.143) is 
— 0.65 /ib. It should be noted that for nuclear targets the 
lowest-lying inelastic channel is the break-up reaction, 
in contrast to photoabsorption on the nucleon where the 
lowest-lying inelastic channel corresponds to single pion 
production. 

No test of the GDH sum rule was hitherto performed 
due to the lack of polarised targets and suitable real- 
photon beams. Only recently, first results of an experi- 
ment on polarised protons in a limited beam energy range 
have been published [Q . Using extrapolations into the un- 
measured regions, Eq. (^) for the proton seems to be sat- 
isfied within the experimental uncertainties. Further real- 
photon experiments are underway at various laboratories 
to extend the energy range of the measurements 

The GDH integral can be generalised to non-zero pho- 
ton virtuality Q 2 in terms of the helicity-dependent virtual 

photon absorption cross sections and 

/*°° A 

Igdh(Q 2 )= [a^(v,Q 2 )-a^(v,Q 2 )} — . (2) 

Jvo v 

The cross section difference appearing in the integrand 
is given by 

Aa = a^-a^ = ^A lFl . (3) 
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In terms of photon-nucleon (nucleus) helicity states this 
relation is valid for any target; in case of the deuteron it 
comprises a mixture of vector and tensor states. Here Ai 
is the photon-nucleon (nucleus) helicity asymmetry, F\ 
the unpolarised nucleon (nucleus) structure function and 
K the virtual-photon flux factor. 

Various generalisations of the GDH integral have been 
considered in the literature. The difference lies in the choice 
made for K. In the notation of Ref. J^] three such gener- 
alisations were considered. In terms of A\ and F\ they 
read: 



Ia(Q 2 



Ib(Q 2 



Ic(Q 2 ) = 



87r 2 a 



87r 2 a 



87r 2 a 



AiFxdx, 



-.A 1 F 1 dx, 



-AiFidx, 



(4) 



(5) 



(0) 



with x — Q 2 /2Mv. Ia corresponds to the case K = v. 
The Gilman notation K = v^J\ + ■y 2 |l7|] has been used 
for 1b while for Ic the Hand convention K = u(l — x) 
]l8| was chosen. They all are numerically close to each 
other in the limits of deep inelastic scattering and real 
photon absorption, but lead to different numerical results 
for intermediate Q 2 . As was pointed out in Ref. J8|, the 
generalisation given in Eq. (ph is most clearly related to 
photoabsorption cross sections. Hence, the generalisation 
used for the figures in this paper is Ib ■ The full numerical 
results will be given for all three prescriptions. 

When considering a nucleon target (spin ^, mass M) 
the photon helicity asymmetry A\ is identical to the longi- 
tudinal virtual-photon asymmetry A\ and the generalised 
GDH integral can be written in terms of the spin structure 
functions g\ and g 2 as: 



GDH 



(Q 2 ) 



8ir 2 a [ Xo gi{x,Q 2 )- 1 2 g 2 {x 1 Q 2 )dx 



M 



K 



(7) 



where g\ and g 2 are the polarised structure functions of 



the nucleon, y 2 = Q 2 /v 2 , xq 



i /2Mv . 



Examining the generalised GDH integral provides a 
way to study the transition from polarised real-photon 
absorption (Q 2 = 0) on the nucleon to polarised deep in- 
elastic lepton scattering (DIS). In other words, it consti- 
tutes an observable that allows the study of the transition 
from the non-perturbative regime at low Q 2 to the per- 
turbative regime at high Q 2 . Since the generalised GDH 
integral is calculated for inelastic reactions, elastic scatter- 
ing is excluded from its calculation. As has been pointed 
out in Ref. || , the elastic contribution to the photon cross 
section becomes the dominant one below Q 2 ~ 0.5 GeV 2 ; 
it has to be taken into account when comparing with twist 
expansions of the first moment of the spin structure func- 
tion g\. In the kinematic region considered in this paper, 
elastic contributions are expected to be small. 

Assuming that the Burkhardt - Cottingham sum rule 



g 2 (x, Q 2 )dx = 



(8) 



holds in good approximation due to the relatively large 
Q 2 values considered in this paper, then Eq. (Q) simplifies 
to 

167T 2 Q 



Ia(Q 



Q 2 



: A(<n 



(9) 



As Q 2 becomes larger, the other generalisations Ib and 
eventually Ic also converge to this value. The first mo- 
ment of the spin structure function gi, J\ = gi{x)dx, is 
predicted to have at large Q 2 only a logarithmic Q 2 depen- 
dence from QCD evolution. Since for the proton rf > 
for higher Q 2 , I GDH must change sign as Q 2 approaches 
zero in order to reach the negative value predicted by the 
GDH sum rule at the real-photon point. The different gen- 
eralisations lead to different values for the expected zero 
crossing needed to connect the negative value predicted by 
Eq. (nj) with the positive value required by measurements 
of -Tfin the DIS limit. For the neutron rf is negative for 
all measured Q 2 . 

The difference of the GDH integral for the proton and 
the neutron, Igdh~^gdh^ IS °f g rea t interest. In the real- 
photon case, the GDH sum rule gives I GDH — Iqdh = 
29/ib, with a sign opposite to what results from mul- 
tipole analyses of meson photoproduction data [^0). In 
the Bjorken limit the difference rf — rf is given by the 
Bjorken sum rule. It can be derived using only current al- 
gebra and isospin symmetry fll|| . This sum rule relates the 
difference of the first moments of g\ and g" at fixed Q 2 to 
the well-measured neutron beta-decay coupling constant 
9a = \g A /gv\ = 1.2670 ±0.0035 @: 



rf-Tf = \-g a -C ns (a s {Q 2 )), 



(10) 



where C ns is the non-singlet QCD correction calculated 
thus far up to O(ck^) in the modified minimal subtrac- 
tion (MS) scheme Experimental verification of the 
Bjorken sum rule at finite Q 2 provides a fundamental test 
of QCD. A measurement of I GDH — Igdh a * large enough 
Q 2 provides such a test. Previous measurements are con- 
sistent with the sum rule when perturbative QCD correc- 
tions are included 

The Q 2 -dependence of the generalised GDH integral 
can be studied separately in the DIS region, characterised 
by large photon-nucleon invariant mass squared W 2 = 
M 2 + 2Mv — Q 2 , and in the nucleon-resonance region 
where W 2 amounts to only a few GeV 2 . Several experi- 
ments measure the generalised GDH integral at low and 
intermediate Q 2 , but cover kinematically only the low- 
W 2 region ]l9[,p0|,pl|. On the other hand, the high-W 2 
contribution to the generalised GDH integral is found to 
be sizeable and essential to any estimate of the total in- 
tegral [p2|,p3"|. Preliminary data from real-photon experi- 
ments at higher energies support this statement ||. The 
kinematics of the HERMES experiment allow the study 
of the Q 2 -development of the generalised GDH integral 
simultaneously in both the nucleon-resonance and DIS re- 
gions. 

In section |^ the experimental setup for data taken 
with a deuteron (proton) target will be described followed 
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by a description of the analysis procedure for both tar- 
gets in section ^. The results for the deuteron nucleus are 
presented in section || together with the proton data re- 
analysed with respect to Ref. |p3[ using an updated value 
for the target polarisation. From these two data sets the 
value of Iq DH is calculated in section [|. Here the assump- 
tion is made that in the kinematical range under consid- 
eration, nuclear effects are small and the deuteron can be 
treated as consisting of two quasi-free nucleons. The re- 
sults on the deuteron nucleus, the proton and the neutron 
are discussed in section ^. From the values on the proton 
and neutron, the proton-neutron difference is calculated 
and compared to the Bjorken sum rule prediction in sec- 
tion ^. A summary of the paper is given in section ||. 



2 Experiment 

HERMES data on the deuteron target were taken in 1998 
to 2000 with a 27.57 GeV beam of longitudinally polarised 
positrons incident on a longitudinally polarised atomic 
Deuterium gas target internal to the HERA storage ring at 
DESY. Data on the proton were taken in 1997 using a lon- 
gitudinally polarised atomic Hydrogen target. The lepton 
beam polarisation was measured continously using Comp- 
ton backscattering of circularly polarised laser light |Q 
|25| . The average beam polarisation for the deuteron (pro- 
ton) data set was 0.55 (0.55) with a fractional systematic 
uncertainty of 2.0% (3.4%). 

The HERMES polarised gas target |2(J consists of po- 
larised atomic D (H) confined in a storage cell. It is fed 
with nuclear-polarised atoms by an atomic-beam source 
based on Stern-Gerlach separation |27j] and provides an 
areal target density of about 2x 10 14 (7x 10 13 ) atoms/cm 2 . 
The nuclear polarisation of atoms and the atomic fraction 
are continously measured with a Breit-Rabi polarime- 
ter Jl^] and a target gas analyser [^| , respectively. The po- 
larisation of the atoms can be flipped within short time in- 
tervals providing both vector-polarisation states and thus 
minimising systematic effects in spin-asymmetry measure- 
ments. The average value of the target polarisation for the 
deuteron (proton) data was 0.85 (0.85) with a fractional 
systematic uncertainty of 3.5 (3.8)%. The value of the pro- 
ton target polarisation used for the data presented in this 
paper has been updated with respect to Ref. making 
use of improved knowledge of samp ling corrections and 
treatment of molecular polarisation [|30f| - The luminosity 
was monitored by detecting Bhabha events using calorime- 
ter detectors close to the beam pipe |3l|. The integrated 
luminosity per nucleon of the deuteron (proton) data set 
was 222 pb- 1 (70 pb" 1 ). 

Scattered positrons, as well as coincident hadrons, were 
detected by the HERMES spectrometer |52). Positrons 
were distinguished from hadrons with an average efficiency 
of 99% and a hadron contamination of less than 1% us- 
ing the information from an electromagnetic calorimeter, 
a transition-radiation detector, a preshower scintillation 
counter and a Cherenkov counter. Only the information 
on the scattered positron was used in this analysis. 



3 Data analysis 

In the following, the analysis procedure used for the deute- 
ron data is given. The analysis procedure and treatment of 
systematic uncertainties have been taken from Refs. |22l 
|33| and are detailed in Ref. [^3|, where the same anal- 
ysis for the proton data was performed. For complete- 
ness, the values and parametrisations used in the latter 
are given below. Note that, compared to Ref. the 
proton data set has been re-analysed in the full kinematic 
range of Ref. [^3| to optimise the binning of the kinemati- 
cally more restricted nucleon-resonance region, where the 
detector acceptance prevents the full coverage over Q 2 . 

The kinematic requirements imposed on the scattered 
positrons in the analysis were identical for both targets. 
The full range in W 2 (1.0 < W 2 < 45 GeV 2 ) was sep- 
arated into nucleon resonance region (1.0 < W 2 < 4.2 
GeV 2 ) and DIS region (4.2 < I^ 2 < 45.0 GeV 2 ). The Q 2 - 
range 1.2 < Q 2 < 12.0 GeV 2 was divided into six bins; 
the same binning as in the proton case was chosen for the 
analysis of the deuteron data and for the subsequent de- 
termination of Iq DH - After applying data quality criteria, 
0.55 (0.13) million events on the deuteron (proton) in the 
nucleon-resonance region and 8.3 (1.4) million events in 
the DIS region were selected. 

For all positrons detected, the angular resolution was 
better than 0.6 mrad, the momentum resolution (aside 
from Brcmsstrahlung tails) better than 1.6% and the Q 2 - 
resolution better than 2.2%. The threshold Cherenkov de- 
tector used in the proton measurement was replaced by 
a Ring-imaging Cherenkov detector |34j] for the data tak- 
ing on the deuteron. The additional amount of material 
led to a slightly worse W 2 -resolution of SW 2 ss 1.0 GeV 2 
for the deuteron as compared to the proton measurement 
(5W 2 w 0.82 GeV 2 ) . Although these ^-resolutions do 
not allow distinguishing individual nucleon resonances, 
the integral measurement in the nucleon-resonance region 
is not degraded. 

The generalised GDH integral Eq. (^) can be re-written 
for any target in terms of the photon-target helicity asym- 
metry A\ and the unpolarised structure function F\. 



87r 2 a f x ° Ax (x, Q 2 )Fi (a;, Q 2 ) dx 



K 



(11) 



where K is the virtual-photon flux factor. 

The cross-section asymmetry A\ for the absorption of 
virtual photons was calculated from the measured cross 
section asymmetry Aii as 



A x 



4l 

D 



rjA 2 



(12) 



For spin-i targets the photon helicity asymmetry Ax is 
identical to the longitudinal virtual photon asymmetry Ax 
and Ai is identical to A 2 . The difference between these two 
asymmetries is relevant for the deuteron target only. Even 
here it is considered to be small in the kinematic region 
examined in this paper and hence will be neglected in the 
following. 



5 



The measured cross section asymmetry Am is given by 



A„ 



(13) 



Here N is the number of detected scattered positrons, L is 
the integrated luminosity corrected for dead time and L p 
is the integrated luminosity corrected for dead time and 
weighted by the product of the beam and target polari- 
sations. The superscript =>(■<=) refers to the orientation 
of the target spin parallel (anti-parallel) to the positron 
beam polarisation. The kinematic factor r\ is given by 



V 



7 (1 - y - tV/4) 

(i-i//2)(i + rV2)' 



(14) 



where y = v/Ef, eam is the inelasticity of the reaction. The 
effective polarisation of the photon D 



D 



y(2-y)(l + 1 2 y/2) 



y 2 (l+ 7 2 )(l-2m 2 /Q 2 ) + 2(l-y 



7 V/4)(l + #) 
(15) 

depends also on R = ulI<^t, the ratio of the absorp- 
tion cross sections for longitudinal and transverse virtual 
photons and the electron mass m e . A 2 is related to lon- 
gitudinal-transverse photon-nucleon interference and is 
not measured in the present experiment. In the DIS re- 
gion Ai can be parametrised in a general form as A2 = 
cMx/ \/Q 2 , where c is a constant determined from a fit 
to the data given in Refs. @,P| as c = °- 20 (°- 53 ) for 
the deuteron (proton) . In the nuclcon- resonance region no 
data are available for the deuteron and A2 — was chosen, 
while for the proton a constant value of A 2 — 0.06 ± 0.16 
was adopted as obtained from SLAC measurements at 
Q 2 = 3 GeV 2 Q. 

Radiative effects for both targets were calculated using 
the codes described in Ref. [[35]. They were found not to 
exceed 7% (4%) of the asymmetry A\ for the deuteron 
(proton). On the integral level they do not exceed 2% and 
were included in the systematic uncertainty. 

The fraction of events smeared from the DIS to the 
nucleon-resonance region and vice versa is evaluated by 
a Monte Carlo simulation of both regions including ra- 
diative and detector effects. Smearing effects in the deep 
inelastic region have been evaluated for all targets follow- 
ing the procedures described in Ref. pj|. The events on 
the deuteron (proton) were simulated using the parametri- 
sation of F2 from Ref. |36| (|37|) for the DIS region, the 
elastic form factors from Ref. jfjq] ( ]39[| ) and the parametri- 
sation of F2 in the nucleon-resonance region from Ref. jlo| 
for both targets. Fig. [l] shows the distribution of experi- 
mental data as a function of W 2 in comparison with the 
simulated events on the deuteron. It is apparent that the 
shape of the simulated distribution agrees well with the 
data. Similar agreement has been found for the proton. 

For the deuteron (proton) case, the relative contami- 
nations from the quasi-elastic (elastic) and deep inelastic 
region in the nucleon-resonance region range from 15% 
(10%) to 3% (2%) and from 11% (7%) to 23% (16%) re- 
spectively, as Q 2 increases from 1.2 GeV 2 to 12.0 GeV 2 . 



a 



5 20000 



515000 



3 

Z 



10000 



5000 




6 2 2 8 
W 2 (GeV ) 



Fig. 1. Comparison of deuteron data with Monte Carlo sim- 
ulations for the nucleon-resonance region as a function of W 2 . 
The total simulated distribution has been normalised to the 
data. The cross-hatched area represents the contribution from 
quasi-elastic scattering, while the lined areas show the contri- 
bution from the nucleon-resonance region (left) and from the 
DIS region (right). The solid line indicates the sum of all sim- 
ulated events and compares favourably with the data points. 
The statistical uncertainties of the data are covered by the 
symbols. 



The fraction of events smeared from the nucleon-resonance 
region to the deep inelastic region ranged from 2.9% (2.5%) 
to 0.5% (0.2%), respectively. Smearing from the elastic 
region to the DIS region can be neglected in the present 
experiment. 

To evaluate the systematic uncertainty from smearing, 
two different assumptions on A\ for the deuteron (pro- 
ton) have been used: a polynomial representation A\ = 
-0.0307 + 0.92a; - 0.28x 2 (power law A x = x 727 ) that 
smoothly extends the DIS behaviour for the asymmetry 
into the nucleon-resonance region |0| ; and for both tar- 
gets a step function (Ai = -0.5 for W 2 < 1.8 GeV 2 and 



At 



1.0 for 1.8 GeV < W < 4.2 GeV ) that is 



suggested by the hypothesis of the possible dominance of 
the P33 -resonance at low W 2 and of the Sn -resonance at 
higher W 2 (see e.g. Ref. p2|). The combined systematic 
uncertainty in the partial integrals from smearing and ra- 
diative effects does not exceed 14% (10%) for the deuteron 
(proton) data. In both cases, smearing gives by far the 
dominant contribution. 
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Fig. 2. The cross section difference Aa obtained for the 
deuteron nucleus as a function of the virtual photon energy 
v for various bins in Q 2 . The dashed curves represent the 
parametrisation used to extrapolate into the unmeasured re- 
gion at high W 2 . 

4 Results for Deuteron and Proton 

The GDH integrals for the deuteron and proton were eval- 
uated following the procedure described in the previous 
section. The nucleon-resonance region and the DIS re- 
gion were treated separately. The large H /2 -range covered 
by the HERMES experiment allows essentially the first 
experimental determination of the complete generalised 
GDH integral for the deuteron, proton and neutron. 

The GDH i nteg ral 1% DH f° r the deuteron was evalu- 
ated using Eq. (O) in both the nucleon-resonance region 
and the DIS region. Here and in the following Iq DH is un- 
derstood as the generalised GDH integral for the deuteron 
nucleus. The unpolarised structure function Ff = Ff (1 + 
7 2 )/ (2x{l + R d )) was calculated in the nucleon-resonance 
region from a modification of the parametrisation of Ff 
given in Ref. ]40|] that accounts for nucleon resonance ex- 
citation assuming R d = <jlI&t to be constant and equal 
to 0.18 in the whole W r2 -range. In the DIS region Ff for 
the deuteron was calculated following a parametrisation 
of F d from Ref. |3(| . In the same kinematic region R was 
chosen according to a fit in Ref. J43| . Note that due to 
cancellations between the R d dependences of Ff and D 
at low y the final result is affected by at most 2% by a 
particular choice of R d . The W 2 dependence of the inte- 
grand Ff / K in the individual bins was fully accounted for 
in the integration. 

The integrand Aa used to calculate Iq DH for the deute- 
ron target is shown in Fig. || as a function of v for the 



Fig. 3. The cross section difference Aa obtained for the proton 
target as a function of the virtual photon energy v for various 
bins in Q 2 . Note that for the upper left plot the first data 
point corresponding to 28.18 ± 6.79 fib at v = 1.8 GeV is off 
scale. The dashed curves represent the parametrisation used to 
extrapolate into the unmeasured region at high W 2 . 

various bins in Q 2 . For the proton case the corresponding 
values for Aa are shown in Fig. ||. In both cases, the ex- 
trapolation into the unmeasured region for W 2 > 45 GeV 2 
was done using a multiple-Reggeon exchange parametri- 
sation |Q for Aa at high energy. The resulting contribu- 
tions are given in Table |lj and range for the deuteron from 
-0.07 lib at Q 2 = 1.5 GeV 2 to 1.53 fib at Q 2 = 6.5 GeV 2 . 
The corresponding contributions for the proton amount 
to about 3.5 \ib for all Q 2 -bins. 

The generalised GDH integrals for the deuteron data, 
calculated in the nucleon-resonance region, in the DIS re- 
gion and over the full W 2 -range, are depicted in Fig. ||. 
The statistical and systematic uncertainties of the full 
Igdh are clearly dominated by the uncertainties in the 
nucleon-resonance region. They are particularly large due 
to the smallness of D and the large size of 77 accentu- 
ating the uncertainties in A d , which amounts to 30% of 
the nucleon-resonance contribution. The systematic un- 
certainty on A d in the DIS region does not contribute 
significantly. The systematic uncertainty for the extrap- 
olation to the unmeasured region at high W 2 of 5% has 
been taken into account. Further sources of systematic 
uncertainties include the beam and target polarisations 
(5.5%), the spectrometer geometry (2.5%), the combined 
smearing and radiative effects (14% of the partial inte- 
grals) and the knowledge of F 2 (5%). The total systematic 
uncertainty of the total GDH integral ranges from 16% at 
Q 2 = 1.5 GeV 2 to 7.5% at Q 2 = 6.5 GeV 2 . For the sys- 
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Fig. 4. The generalised GDH integral Iq DH for the deuteron 
nucleus, shown as a function of Q 2 for the three kinematic 
regions considered: nucleon-resonance region (triangles), DIS 
region (squares), and full VK 2 -region (circles) including extrap- 
olation to the unmeasured part. The error bars show the statis- 
tical uncertainties. The solid curve is taken from Ref. and 
represents a prediction for the full Iq DH . The dashed curve rep- 
resents a model for the nucleon-resonance region from Ref. Ma] . 
The systematic uncertainties of the full integral are given as a 
band; the hatched area inside represents the systematic uncer- 
tainty of the nucleon-resonance region alone. Note that some 
data points are slightly shifted for better visibility. 



tematic uncertainties of the nucleon-resonance and DIS 
regions, independent sources of systematic uncertainties 
were added in quadrature, while the systematic uncertain- 
ties stemming from smearing effects and the knowledge of 
were added linearly. Only smearing effects from the 
quasi-elastic region to the measured range of W 2 had to 
be taken into account for the total integral, thus reduc- 
ing considerably its systematic uncertainty due to smear- 
ing compared to the integrals calculated separately in the 
nucleon-resonance and DIS regions. 

The generalised GDH integrals for the proton data re- 
analysed using an updated value of the target polarisa- 
tion, calculated in the nucleon-resonance region, in the 
DIS region and over the full IF 2 -range, are depicted in 
Fig. [| For both targets, the contribution of the nucleon- 
resonance region decreases faster than that of the DIS 



region as a function of Q . The latter dominates /; 
and I P GDH for Q 2 > 3.0 
at the lowest measured 
tribution shown in Fig. 



GDH 

GeV 2 and remains sizeable even 
Q 2 . The nucleon-resonance con- 
and Fig. || respectively for the 



deuteron and proton is compared to a curve representing 



Fig. 5. The generalised GDH integral Iq DH for the proton us- 
ing the same notation as in Fig. ^. The solid curve is taken from 
Ref. jlfi), the dashed curve follows the model of Ref. Q. The 
data were published earlier in Ref. [S] , but are re-analysed for 
the present paper using improved knowledge on the target po- 
larisation. The band representing the systematic uncertainties 
is given for the convention defined in Fig. [l|. 

the prediction of the model of Ref. This model is 

based on a helicity-dependent sum over the first, second 
and third nucleon-resonance regions using experimental 
resonance parameters, but assuming infinitely narrow res- 
onances. The threshold region was taken into account in 
the first nucleon-resonance region. Within the total ex- 
perimental uncertainties this model describes the data. 

The data for the full integral on both the deuteron 
and the proton target are compared to a model based 
on the leading twist <5 2 -evolution of the first moments 
of the two polarised structure functions g\ and gi with- 
out consideration of any explicit nucleon-resonance con- 
tribution p6,47|. In this model, the low-Q 2 behaviour of 
32 is governed by the Q 2 -dependence of a linear combi- 
nation of the electric and magnetic Sachs form factors. 
The model predicts the shape. It predictions thus depend 
on the experimental value for Z~i at asymptotically large 
Q 2 and on the Q 2 -dependence of the Sachs form factors 
at low Q 2 . The normalisations for asymptotically large 
Q 2 was taken from the present data. They were evalu- 
ated from the Q 2 -dependencies of Iqd H with the l/Q 2 - 
dependence expected from leading twist devided out. Fit- 
ted by straight lines the results are rf — 0.129 ±0.006 and 
If = -0.030 ± 0.007 (cf. Fig. |) . The parameterisation 
of the form factors was taken from Ref. |4^] . The model 
describes the data on the proton very well. No explicit 
prediction for the deuteron is given; thus the deuteron is 



8 




6 2 2 8 

Q 2 (GeV 2 ) 



Fig. 6. The generalised GDH integral Iqdh f° r the neutron 
obtained from the deuteron and proton using the same notation 
as in Fig. ^] for the symbols and theoretical curves. 



modelled as the sum of proton and neutron. Nevertheless, 
the model prediction also agrees well the deuteron data. 



5 Neutron results from Deuteron and Proton 

The extraction of the generalised GDH integral for the 
neutron from data taken on the deuteron and the proton 
requires nuclear effects such as Fermi motion and the de- 
polarising effect of the D-state to be taken into account. 
These questions were addressed in Ref. |48 . Following 
their model, the integral Iq D h for the neutron was calcu- 
lated from the results Iqbh 011 the deuteron, as obtained 
in this analysis, and those on the proton Iqxih re-analysed 
following the procedure detailed in Ref. [ f23[ : 



jn 

1 GDH 



Td 

1 GDH 

1 - 1.5w d 



I p 



GDH' 



(16) 



Here u d = 0.050 ± 0.010 @ is the probability of the 
deuteron to be in a D-state. It has been shown in Ref. p8| 
that although the uncertainties in the structure functions 
in the integrand of Eq. ( |li"| ) may be large, the resulting 
contribution to the systematic uncertainty for the inte- 
gral Iq DH due to nuclear effects does not exceed 3%. This 
combined with the uncertainty in Ud leads to an additional 
systematic uncertainty on Iq DH of 4%. No further as- 
sumptions, in particular not on i 7 ^, A% and R, are needed 
to derive the generalised GDH integral for the neutron 
using Eq. (HO) . 



In the real-photon case, the application of Eq. (|16|) is 
not straightforward, since significant contributions from 
photodisintegration and coherent photoproduction must 
be taken into account Q . For virtual photons three dif- 
ferent regions can be distinguished. For Q 2 > 1 GeV 2 the 
generalised GDH integral can be described by the spin 
structure functions g\ and g2 taking twist-2 and twist- 
3 contributions into account (cf. Eq. (0)). Higher twist 
contributions are suppressed by powers of l/Q 2 ; Eq. (|l^ ) 
holds. For Ql < Q 2 < 1 GeV 2 , Eq. (Jl|) holds, but the 
GDH-integrals for proton and neutron deviate substan- 
tially from the g\ and g2 (twist-2 and twist-3) contribu- 
tions. In this region, the Operator Product Expansion has 
already broken down and one has to resort to model as- 
sumptions like those of Refs. |4(|[47]]. The relevant scale 
Qq was estimated in Ref. [[[l| to be Qq ~ to 2 . Finally, for 
Q 2 < Qq Eq. (fl6| ) is no longer valid. 

In Fig. fj the results for I GDH obtained from the deu- 
teron and proton data in three VF 2 -regions are shown 
together with model predictions following Ref. |Q and 
Ref. f46| . As in the proton case, the contribution from the 
nucleon resonance region decreases faster with increasing 
Q 2 compared to the contribution from the DIS region. 
The data are well described by the resonance model. The 
contribution from the extrapolation to high W 2 is dom- 
inant for Q 2 > 2.0 GeV 2 and remains sizeable down to 
the lowest measured Q 2 (cf. Table [l]). In agreement with 
measurements of the neutron spin structure function g" 
and as expected from recent measurements of polarised 
quark distributions |p2| , Iq DH is negative and of smaller 
absolute size than the proton value. Within the total ex- 
perimental uncertainties, the model prediction of Ref. [fl6) 
agrees well with the neutron data for the full integral. 

Results on Iq DH were also obtained from a previous 
measurement on a 3 He target p2| . The neutron asymme- 
try was obtained from the 3 He asymmetry taking into ac- 
count nuclear effects, the relative polarisation of the neu- 
tron and two protons, as well as a fit to the data for Af[. 
Note that the lower 14 72 -limit for the data taken on 3 He 
was 4 GeV 2 and thus slightly different from the cut at 4.2 
GeV 2 used in the deuteron analysis. Both data sets are 
shown in Fig. ^ and agree within their respective uncer- 
tainties. 



6 Discussion of results 

In Table || the final results are presented for the full gen- 
eralised GDH integrals on the deuteron, the proton and 
the neutron in bins of Q 2 and for the three generalisations 
(Eqs. |], ||, H) considered in the literature. Significant dif- 
ferences between the integral values in various generalisa- 
tions are observed. While I a and Ib remain comparable 
at the measured Q 2 values due to the smallness of 7 2 , 
the weighting introduced by the Hand notation in Iq 
leads to sizeable differences. The results will be discussed 
referring to the generalisation Ib ■ 

As mentioned above, the final results for the proton 
and for the neutron, as presented in Fig. || and Fig. |(], 
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Fig. 7. The generalised GDH integral Iq DH in the DIS region 
for the neutron (squares) obtained combining deuteron and 
proton data, shown in comparison to the results obtained on 
3 He (triangles) |22| . The lower WMimit for the latter was 



W > 4.0 GeV 2 while for the former W > 4.2 GeV 2 was 
used. The error bars represent the statistical uncertainties. The 
systematic uncertainties for each data set are given as error 
bands (deuteron top, 3 He bottom). 



show that the contribution of the nucleon resonance region 
to the full generalised GDH integral is small for Q 2 > 3 
GeV 2 and the contribution from the DIS region remains 
sizeable down to the lowest measured Q 2 values. Numeri- 
cal values following the generalisation Ig are given in Ta- 
ble 0. 

The results for the full generalised GDH integrals on 
the deuteron, proton and neutron are shown together in 
Fig. ||. As noted above, they agree well within the total 
uncertainties with a prediction of Ref. based on the 
leading twist Q 2 -evolution of the two polarised structure 
functions g\ and gi without consideration of any explicit 
nucleon-resonance contribution. Although in the neutron 
case the poorer knowledge of input data for this model 
leads to a larger uncertainty, a similar description of the 
data compared to the proton case is achieved. No turn- 
over is observed in the measured range that would be re- 
quired for the generalised GDH integral on the proton or 
the deuteron to meet the GDH sum rule predictions at 
Q 2 = 0. Preliminary data from Ref. |L9],|2(| and a recent 
theoretical evaluation indicate that this sign change hap- 
pens at a value of Q 2 much lower than the range consid- 
ered in this analysis J53fl . 

At large Q 2 the generalised GDH integral is connected 
to the first moment of the spin structure function g\ (Eq. 



Fig. 8. The Q 2 -dependence of the generalised GDH integrals, 
calculated over the full W 2 -region, for the deuteron nucleus 
(squares), proton (circles) and neutron (triangles). The latter 
was obtained from the deuteron and proton data. The curves 
shown are the predictions for the various targets according to 
Ref. p6| ]. The error bars represent the statistical uncertain- 
ties. The bands represent the systematic uncertainties (open: 
neutron, lined: deuteron, cross-hatched: proton). The open 
symbols at Q 2 — 5 GeV 2 represent the measurements from 
Ref. fjj (shifted to the left) and Ref. (shifted to the right) 
on proton and neutron. The stars represent the three highest 
Q 2 bins of the neutron measurement from Ref. including 
an extrapolation for the unmeasured DIS region. 



@). For Q 2 > 3 GeV 2 the generalised GDH integral is 
completely dominated by the DIS region. The data pre- 
sented in this paper agree with the most recent values for 
the first moments of the spin structure functions measured 
on the proton by E-155 (E-143) rf = 0.118 ± 0.004 ± 
0.007 (rf = 0.129 ± 0.003 ± 0.010) and on the neu- 
tron measured by E-155 (E-143) If 1 = -0.058 ± 0.005 ± 
0.008 (r™ = -0.034 ± 0.007 ± 0.016) evaluated at Q 2 = 
5 GeV 2 0,0. These values correspond to I GD h(Q 2 = 
5 GeV 2 ) = 10.59 ± 0.36 ± 0.73 (11.85 ± 0.28 ± 0.92) fib 
and Igdh(Q 2 = 5 Gc V 2 ) = -5.21±0.45±0.72 (-3.12± 
0.64 ± 1.47) for the proton and the neutron respec- 
tively, and are shown together with the present data in 
Fig. | 

The Q 2 -behaviour of Igdh for proton and neutron can 
be more clearly studied when dividing out from Igdh the 
l/Q 2 -dependence expected from leading twist. According 
to Eq. (||), Igdh is then expected to show a logarith- 
mic Q 2 -dependence, similar to i^. The result is shown 
in Fig. |p. Any contributions from resonance form factors 
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Fig. 9. The Q 2 -dependence of the generalised GDH integrals 
for the proton (filled circles) and neutron(open circles) after the 
leading-twist dependence, Q 2 / (I6iv 2 a) , has been divided out. 
The error bars represent the statistical errors. The systematic 
uncertainties are represented by the respective error bands. 
The dash-dotted lines are straight line fits to the data. 



or higher-twist contributions should become visible as a 
deviation from a flat line. Considering the statistical and 
systematic uncertainties of the present measurement, no 
deviation from a leading-twist behaviour can be seen. In 
other words, the experimental data obtained for the pro- 
ton and the neutron are consistent with the naive expec- 
tation that the l/Q 2 expansion is a good approximation 
down to Q 2 ~ 2 GeV 2 . As discussed in the introduction 
the elastic part (x = 1) has to be included for a complete 
comparison to a twist expansion of r± . However, this is not 
a relevant contribution in the kinematic range considered. 



7 Combined results for the Proton and the 
Neutron 

The data obtained for the proton and neutron over a large 
range in Q 2 and W 2 offer a unique possibility to evaluate 
the proton-neutron difference of the generalised GDH in- 
tegral. This difference is shown in Fig. [l(]. It is expected 
to be less sensitive to higher twist effects or contributions 
from nucleon resonances. Within the measured Q 2 -range 
no turn-over at low Q 2 required to meet the GDH sum 
rule prediction of 29 fib for Q 2 = is observed. The data 
fall off as l/Q 2 indicating that leading twist dominates 
as expected. A fit to the data using c/Q 2 where c is a 
constant is shown together with the data in Fig. [l(]. For 



Fig. 10. The Q 2 -dependence of the generalised GDH integral 
for the proton-neutron difference. The dotted curve represent 
the prediction from Ref. ^7|. The dashed curve represents a 
simple l/Q 2 fit to the data leading to I^ H {Q 2 = 5 GeV 2 ) of 
14.3 ±0.9 ± 1.3 fib. For large Q 2 this difference is expected to 
obey the Bjorken sum rule. The open symbols represent the 
measurements of the Bjorken sum rule from Refs. [jl4J (square) 
and jlH) (circle). They are slightly shifted in Q 2 for clearer 
representation. At Q 2 — the GDH sum rule gives a value 
of 29 fib. The error bars represent the combined statistical 
uncertainties; the band indicates the systematic uncertainty. 



the constant c/(167r 2 a) = 0.159 ± 0.009 is found, leading 



jn 
In 



GDH 



14.3 ± 0.9 ± 1.3 fib at Q 2 



to a value I GDH 

5 GeV 2 . This result is in agreement with an experimental 
determination of the Bjorken sum rule by E-155 (E-143), 
If - r{ 1 = 0.176 ± 0.003 ± 0.007 (0.164 ± 0.008 ± 0.020) 
[ ^4|[l5[ | within the respective experimental uncertainties. 
At Q 2 — 5 GeV 2 these values correspond to I GDH — 

J GDH = 15 - 76 ± °- 27 ± 0- 63 V b ( 14 - 95 ± °- 72 ± 1-79 fib). 
Within errors, the value for the proton-neutron difference 
measured in this analysis is also in agreement with the 
Bjorken-sum-rule prediction of 0.182 ± 0.005 correspond- 
ing to I P GDH - I GDH = 16.33 ± 0.45 fib. The Bjorken 
sum rule was evaluated using Eq. (10) at Q 2 — 5 GeV 2 
including corrections up to third-order in a s and with 
a s (M z ) = 0.118. 



8 Summary 

The generalised GDH integral has been determined for a 
deuteron target in the kinematic region 1.2 < Q 2 < 12.0 
GeV 2 and 1 < W 2 < 45 GeV 2 in this paper. Using an 
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Table 1. The generalised GDH integral Ib for the deuteron, proton and neutron given in [ib per nucleus for the full W /2 -region, 
the nucleon resonance region and the DIS region and various values of Q 2 in GeV 2 . The statistical and systematic uncertainties 
are given. The last three rows show the integral values for the extrapolation to higher W 2 , for the deuteron, proton and neutron. 
Note that the systematic error on the total integral is reduced compared to the integral in the nucleon-resonance and DIS regions, 
since smearing between these regions is not to be taken into account. 



Q 2 


1.5 


2.1 


2.7 


3.5 


4.5 


6.5 




rd 

Hot 
jd 

-*res 
rd 

J DIS 


24.2±3.9 ± 3.9 
11.4±3.6 ±4.0 
12.9±1.3 ± 1.3 


21.9±2.0±3.2 
9.3±1.8±3.3 
12.3±1.0± 1.3 


15.5±1.3 ± 2.4 
4.2±1.1 ± 2.2 
10.6±0.8± 1.3 


12.6±0.7± 1.7 
3.4±0.5 ± 1.6 
8.3±0.5 ± 0.9 


9.0±0.5± 1.2 
0.9±0.3±0.9 
6.9±0.4±0.7 


6.4±0.3±0.5 
0.1±0.1 ± 0.1 
4.8±0.3 ± 0.5 




-'tot 

TV 
-*res 

TV 

I DIS 


42.7±5.7±3.8 
22.2±5.4±4.2 
17.2±1.9 ± 1.1 


28.8±2.9±2.1 
10.7±2.6±2.1 
14.6±1.3± 1.0 


21.7±1.9± 1.2 
5.0±1.5±0.9 
13.2±1.0 ± 1.0 


17.1±1.0±0.9 
2.5±0.7±0.4 
11.1±0.7±0.8 


12.7±0.6 ± 0.6 
1.3±0.4±0.2 
7.9±0.5±0.6 


8.9±0.3±0.5 
0.1±0.1 ± 0.1 
5.5±0.3±0.5 



rn 

J-tot 

rn 
^res 

rn 

J DIS 


-16.5±7.1 ±5.7 
-9.8±6.7±6.1 
-3.3±2.3± 1.8 


-5.1±3.6±4.0 
-0.6±3.2±4.0 
-1.3±1.7± 1.7 


-5.0±2.4±2.9 
-0.6±1.9 ± 2.6 
-1.6±1.4± 1.7 


-3.5±1.3±2.0 
1.2±0.9 ± 1.7 
-2.1±0.9± 1.3 


-3.0±0.8± 1.4 
-0.3±0.5± 1.0 
-0.4±0.6± 1.0 


-1.9±0.4±0.7 
0.1±0.1 ±0.2 
-0.2±0.4±0.7 


jd 

1 unm 


-0.1 


0.3 


0.7 


0.9 


1.2 


1.5 


JV 


3.3 


3.5 


3.5 


3.5 


3.5 


3.3 


jn 
unm 


-3.4 


-3.2 


-2.8 


-2.6 


-2.3 


-1.8 



Table 2. The generalised GDH integral for the deuteron, proton and neutron, calculated for the full W /2 -range, given in fib 
per nucleus for the generalisations given in Eqs. ^, ^, [] following Ref. ji| and various values of Q 2 in GeV 2 . The statistical and 
systematic uncertainties are given. 



Q 2 




1.5 


2.1 


2.7 


3.5 


4.5 


6.5 




rd 
1 A 

rd 
L B 

rd 
J C 


27 .4 
24.2 
44.9 


±4.7 ±4.5 
±3.9 ±3.9 
±9.3 ±7.3 


24.5 ±2.4 ± 3.5 
21.9 ±2.0 ± 3.2 
43.3 ±5.2 ± 6.3 


16.7 ±1.6 ±2.6 
15.5 ±1.3 ±2.4 
28.3 ±3.6 ±4.4 


13.5 ±0.8 ± 1.8 

12.6 ±0.7 ± 1.7 
24.5 ±1.9 ±3.3 


9.4 ±0.5 ±1.2 
9.0 ±0.5 ± 1.2 
14.9 ±1.3 ±2.0 


6.6 ±0.3 ±0.5 
6.4 ±0.3 ±0.5 
9.4 ±0.4 ±0.8 




T P 
A 
TV 
1 B 

T p 
1 C 


48.7 
42.7 
81.7 


±6.9 ±4.3 
±5.7 ±3.8 
±14.0 ±7.4 


31.7 ±3.4 ±2.3 

28.8 ±2.9 ± 2.1 
53.4 ±7.7 ±3.9 


23.2 ±2.2 ± 1.3 
21.7 ±1.9 ±1.2 
37.4 ±5.3 ±2.1 


17.9 ±1.1 ± 1.0 
17.1 ±1.0 ±0.9 
27.8 ±2.8 ± 1.6 


13.3 ±0.7 ±0.6 
12.7 ±0.6 ±0.6 

20.4 ±1.9 ±0.9 


9.0 ±0.3 ±0.5 
8.9 ±0.3 ±0.5 
11.8 ±0.4 ±0.7 




rn 

*A 
rn 
IB 
rn 
1 C 


-19.0 ±8.6 ±6.5 
-16.5 ±7.1 ±5.7 
-33.1 ±17.2 ± 10.8 


-5.2 ±4.3 ±4.5 
-5.1 ±3.6 ±4.0 
-6.6 ±9.5 ±7.8 


-5.1 ±2.8 ±3.1 
-5.0 ±2.3 ±2.9 
-6.6 ±6.6 ±5.1 


-3.3 ±1.4 ±2.2 
-3.5 ±1.3 ± 2.0 
-1.4 ±3.5 ±3.8 


-3.1 ±0.9 ± 1.5 
-3.0 ±0.8 ± 1.4 
-4.3 ±2.3 ± 2.3 


-1.9 ±0.4 ±0.7 
-2.0 ±0.4 ±0.7 
-1.6 ±0.6 ± 1.0 



updated value of the target polarisation, a correspond- 
ing re-analysis of the HERMES proton data as compared 
to Ref. |2j| was performed. In both cases the I4 y2 -range 
was separated at W 2 — 4.2 GeV 2 into a region where 
the nucleon resonances dominate and into the DIS region. 
Combining both data sets, the generalised GDH integral 
for the neutron was calculated in the same kinematic re- 
gions. These neutron results obtained from the deuteron 
agree with those obtained earlier on a 3 He target in the 
same kinematic region. 

Altogether, a complete set of measurements of the gen- 
eralised GDH integrals for the deuteron, proton and neu- 
tron is available. In all three cases the nucleon-resonance 
contribution to the generalised integral decreases rapidly 
with increasing Q 2 and the contribution from the DIS re- 
gion is still sizeable even at the lowest measured Q 2 , em- 
phasising the importance of measuring the GDH integral 



over a large W 2 -r&nge. At larger Q 2 the measured values 
agree well with measurements of the first moments of the 
spin structure function g\. 

For the generalised GDH integrals of the proton, the 
neutron and the proton-neutron difference, the Q 2 -de- 
pendence is in agreement with a leading-twist behaviour; 
within the experimental uncertainties it exhibits no signif- 
icant contribution from either higher twist or resonance 
form factors. The proton-neutron difference is in agree- 
ment with the Bj or ken-sum-rule prediction evaluated at 
Q 2 = 5 GeV 2 within the experimental uncertainties. 
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